Effective Temperature of Red Blood Cell Membrane Fluctuations 
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Biologically driven non-equilibrium fluctuations are often characterized by their non-Gaussianity 
or by an "effective temperature", which is frequency dependent and higher than the ambient tem- 
perature. We address these two measures theoretically by examining a randomly kicked "particle" , 
with a variable number of kicking "motors" , and show how these two indicators of non-equilibrium 
behavior can contradict. Our results are compared with new experiments on shape fluctuations of 
red-blood cell membranes, and demonstrate how the physical nature of the motors in this system 
can be revealed using these global measures of non-equilibrium. 
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Experimental and theoretical studies of biological sys- 
tems confront the issue of active elements which give rise 
to fluctuations that distinguish living matter from inan- 
imate soft-matter systems. Examples range from molec- 
ular motors in the cytoskeleton [l| and active membrane 
pumps [2[ to larger scale objects such as swimming bac- 
teria Q . As with other non-equilibrium systems such as 
driven granular matter, it is unclear how to define useful 
measures for non-equilibrium "activity" . Spontaneous 
fluctuations may be compared with the response of the 
system to small external perturbations, to define an ef- 
fective temperature T e g using the fluctuation-dissipation 
(FD) formalism In most cases T c g is frequency de- 
pendent (unlike the thermal case) , and is larger than the 
ambient temperature. These features quantify the non- 
thermal activity in the system. Another parameter that 
is useful for characterizing deviations from equilibrium is 
the non-Gaussianity k of the distribution function 

In biological systems the nature of the microscopic ac- 
tive elements is difficult to study directly. We demon- 
strate how global statistical measures of activity can be 
used to extract qualitative and quantitative properties of 
the underlying molecular motors. We are motivated by 
living cells, in which the activity is induced by multiple 
motors throughout the system that are directly coupled 
to local degrees of freedom. It is important to explore the 
effects of the number of motors and their level of activity 
on the non-equilibrium nature of the fluctuations Q . 

One well-studied active system is the membrane of the 
red-blood cell (RBC). However, the nature of the molec- 
ular motors in this system is still far from being well 



understood. We question how T e g and k correlate (or 
not) with each other, and which properties of the non- 
equilibrium system do they probe. We focus on RBC 
and present new experimental measurements, but obtain 
results on the non-equilibrium statistical mechanics of 
active systems in general. We introduce a simple model 
of a randomly kicked "particle" , with a variable num- 
ber of kicking "motors" (force producing elements). Our 
generalized particle and motors may represent different 
objects in different systems, and we will be more spe- 
cific when comparing to RBC experiments. We compute 
both T e ff and k, and identify situations in which these 
two non-equilibrium indicators contradict. 

Model. We consider the following overdamped 
Langevin equation for the velocity v = x, 



v = -Xv + f T + Ja + Jr- 



(1) 
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A is the damping coefficient. The thermal force fr(t) is 
an uncorrelated Gaussian white noise: (/t(£)/t(^')) = 
2XTBS(t — t'), with Tb the ambient temperature, and 
Boltzmann's constant set to ks = 1. For the active force 
fA{t) we assume that each of the N m motors produces 
pulses of force ±/o, of duration At. We assume sym- 
metry with respect to the force direction, which is mo- 
tivated by experimental observations of nearly symmet- 
ric active fluctuations of cells @ . While the pulses turn 
on randomly as a Poisson process with an average wait- 
ing time r, unless otherwise stated, we take a constant 
pulse length At. The power stroke of molecular motors 
is a realization for such a relatively well-defined impulse 
length [I0| • We will also consider stochastic pulse lengths 
with an arbitrary distribution P(Ar), and show that if 
-P(At) is Poissonian, the force correlations reduce to the 
shot- noise form studied in [ll| ■ 

To measure the linear response Xxx of the particle po- 
sition x, we apply a small force fn = Foe lult , and find 
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that the position is perturbed as (Sx(t)) = Xxxi^Foe 1 "* , 
with Xxx(u) = (uj(i\ — u)) , Sx being the small change 
in position, irrespective of the active force [l2| ■ 

Due to the linearity of Eq. (UJ), and since fx and Ja 
are uncorrelated, the velocity autocorrelation is [l2| 



S vv (uj) 



2XT E 
A 2 + u 



2N m v 2 X 2 [1 - cos(wAr)] 
(t + At)u 2 (A 2 +w 2 ) 



(2) 



where Vq = fo/X is the asymptotic velocity that the par- 
ticle approaches due to the activity of a single motor. 

Effective Temperature. In equilibrium, the FD theorem 
connects the imaginary part of \xx to the position au- 
tocorrelation, S xx (lo) = -U)~ 2 S vv (u), by: lia[xxx{^)] = 
Tpj;S xx (uj). For non-equilibrium steady states we define a 
frequency-dependent effective temperature 



}S xx (cj) 



21m \Xxx(w)] 



= T E 



N m v 2 X [1 - cos(wAt)] 
(t + At)lu 2 



.(3) 



Note that the time between pulses enters only through 
the density of pulses per unit time. Hence this result 
does not depend on the distribution of waiting times be- 
tween pulses, but only on its average, r. Moreover, our 
results may be extended to a stochastic pulse length 12 1 . 
In particular, for Poissonian At we obtain shot-noise 
force correlations (/a(*)/a(0)) = (At)" 1 cxp (-t/(Ar)) 
andT cff = T B + 7V mV 2 A(T + (Ar))- 1 (u 2 + (At)- 2 )' 1 /2. 
Note that alternative definitions of the effective temper- 
ature have appeared in the context of granular gases [l3| . 
Unlike those systems, here we do not identify a non- 
equilibrium situation where T c ff (w) = const. In the high- 
frequency limit the active contribution vanishes, so that 
T c ff — > Tg. Around lo = 1/Ar, T e g rises and as u — > 0, 
it approaches a constant value in the low frequency limit 
(w « At" 1 ): T cff (0) = Tb + N m v 2 XAT 2 /{r + At). 

Velocity Distribution and Kurtosis. In Fig. [IJa, we 
plot the velocity distribution function P(v) from nu- 
merical simulations. The distribution is highly non- 
Gaussian for a single motor. Interestingly, for small 
AAt and in the presence of multiple motors, P{v) 
may retain its Gaussian form, even though the sys- 
tem is very far from equilibrium, as can be seen in 
the strong frequency dependence of T c ff(cj), and in the 



fact that 



is significantly larger than Tb- For our 



model we can exactly calculate [12(: (v 2 ) = Tb + 
[N m v 2 (AAt + e- AAr - l)] / [A (r + At)]. Let us em- 
phasize that (v 2 ) ^ T c ff(0), but rather in the limit 
AAt -> we find that (v 2 ) = 2T cS (0). 

We measure the non-Gaussianity by the kurtosis, K = 
(v 4 ) I (v 2 ) 2 , which we plot in Fig. [TJd as a function of the 
number of motors and the activity quantified by the prob- 
ability of a single motor to be on: p on = At/(t + At). 
To calculate (v 4 ) for a single motor, as long as At ^ 1, 
we ignore overlaps between the contributions of consec- 
utive pulses (l2j . A simple model which works rather 
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FIG. 1: (a) Velocity distribution function. Thermal (dash- 
dot straight line with slope T _1 visible in the inset), N m = 1 
(solid) and 10 (dashed). A = 50, At = 0.1, r = 0.15, T B = 



10" 



«o = 0.02. Note the peak at vo for N m = 1. Dotted 



line: N m = 10 using T B = 10~ 7 and A = 150. The kurtosis 
for the active cases are: k/kg = 0.85,0.99,0.96, respectively, 
(b) Kurtosis vs motor activity (varied by changing t) from 
numerical simulations (symbols), compared to the analytic 
expression ignoring pulse overlaps (dotted lines), and to the 
model of shifted Gaussians (solid lines). 



well at all number of motors, approximates the distri- 
bution as a sum of shifted thermal Gaussians |l2|. As 
long as AAt 3> 1 , the contribution to the velocity distri- 
bution due to the rise and decay before and after each 
pulse is small, and this model gives a very good descrip- 
tion (see Fig. [T|b). In the opposite limit of At <C 1 and 
AAt <C 1 the velocity distribution approaches a Gaus- 
sian. The value of the kurtosis measures the spread of 
the distribution; larger values correspond to a distribu- 
tion that is wider than a Gaussian, and vice versa. 

The most outstanding result is the non-monotonic de- 
pendence of k on p on . Compared to a dead system 
(Pon — 0), as the motor activity is turned on, and when 
the active velocity Vq > Tb, the velocity distribution gets 
more populated at higher values, hence n increases. In 
the other limit of p on — > 1, k is necessarily smaller than 
the Gaussian value kg = 3, since it is a contribution of 
shifted Gaussians. From these two limits we conclude 
that k is necessarily a non-monotonic function of p on . In 
fact, k can retain its kg value even for a non-equilibrium 
system (p on > 0). For the distributions shown in Fig.QJi, 
k is close to kg, except for the single motor, even when 
the distribution is visibly non-Gaussian (N m = 10). As 
a function of N mi at small p on the deviation from kg 
increases with N m (Fig. [TJd), since the high velocity tails 
are more populated. In the other limit of p on — > 1, the 
distribution approaches a Gaussian with increasing N m , 
which is a manifestation of the central limit theorem. 

Comparing with T e g we find that both measures of 
non-equilibrium behavior increase with increasing activ- 
ity in the p on — > limit, while as p on —> 1 they contradict. 
Both N m and p on are simply multiplied to give the ampli- 
tude of the active contribution to T e g, while k is a more 
complicated function of these two parameters. 

Experiments. The activity of the RBC membrane was 
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FIG. 2: (a) Relative kurtosis vs wavevector: Squares (circles) 
are for natural (starved for 6hrs) RBC. k,g(<i) was extracted 
from the data on ATP-depleted cells, starved for 24hrs. (b) 
Calculated dependence of k from our model, mapped to q- 
space (using p on = 0.07). The dimensionless q is determined 
by varying the number of motors {N m = 1, 2. ..30) 
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FIG. 3: Effective temperature for a free membrane of lateral 
size L = 8^im, driven by the direct force (solid line) and 
curvature force (dashed line). Inset: L = 50/xm. 



recently measured in two different experiments, which 
found indications for non-equilibrium fluctuations when 
the chemical energy source of ATP is available. Be- 
fore comparing with our model we note that the mem- 
brane undulations may be described by the following 
over-damped analogue of Eq. ([!} [l6| , 



-Xqh q 



O q [F T {q,t) + F A {q,t)} 



(4) 



where h q is the amplitude of the membrane deflection 
at wavevector q, X q = O q (Rq i + oq 2 ) is the response of 
the membrane due to the elastic restoring forces of cur- 
vature and tension (with bending modulus R and mem- 
brane tension a), O q = (Arjq)^ 1 is the Oseen interaction 
kernel for a flat membrane in free fluid and rj the viscos- 
ity of the surrounding fluid. The thermal force satisfies 
(F T (q,t)F T (-q,t')) = 2T B Q- 1 5(t - t'), and F A (q,t) is 
the Fourier transform of the active force. For the ac- 
tive forces we consider two direct force and a 
curvature- force both with shot-noise correlations. 

The first experiment 0] measured the spatial depen- 
dence of the membrane fluctuations, and extracted the 
probability distribution P(h q ), from which the kurto- 
sis was obtained. Here, k > kq was found for ATP- 
containing cells. In Fig. [2£i we present new data showing 
that k increases with q and with the ATP concentra- 
tion [r3 |. Comparing these observations with our model 
(Fig.[T]) this indicates that the RBC has p on close to zero, 
which means that r 3> At. 

Next, we compute the ^-dependence of the kurtosis us- 
ing our single-particle model. We map each mode q of 
the membrane to a single particle as follows; the number 
of motors that act on the membrane area involved in the 
motion of mode q is given by: N m oc q~ 2 (the number 
of motors in the membrane area of wavelength 2tt /q, as- 
suming they are uniformly distributed on the membrane), 
Fa oc q°,q 2 (direct and curvature force respectively 
and (Fj) oc q [l2|. Figure EJd shows that our calculation 
predicts that k — > kq as q — > 0, in agreement with the 
experiment (Fig. (5^.). For small q, N m increases, and we 



are in the low p on regime. Note that if the RBC had 
large p on , k would decrease with increasing q. The peak 
in the experimental data may indicate the wavelength 
corresponding to a single active unit ("motor") in the 
RBC cytoskeleton . Note that using our single-particle 
model for the dynamics of an extended object such as the 
membrane is only a qualitative approximation. 

Thus, the motors are distributed throughout the RBC 
membrane, and each motor has a long recovery time. 
These findings agree with the proposed mechanism of 
membrane fluctuations jTijl; ATP induces the release of 
membrane-anchored filaments, the release event (At) is 
fast compared to the time it takes the released polymer 
to find its anchor on the membrane and re-attach (r). 

Another experiment 15[ measured the frequency de- 
pendence of the height fluctuations at a single point on 
the RBC membrane, and found a 3-7 fold increase in 
low frequency (/ < 10Hz) fluctuations compared to cells 
depleted of ATP. The way to decouple the ATP-induced 
changes to the elastic moduli [14( from the increase in T e s 
is to measure the response in addition to the fluctuations, 
and this awaits future experiments. In Fig. [3] we plot the 
calculated effective temperature of the system, as defined 
by Eq. (j3]) 12j. T e g approaches Tb for large frequencies 
uj ^> t _1 , and increases for frequencies uj < At -1 . There 
is even a peak in T e g for the curvature-force. The val- 
ues of T e ff reach up to 10Tb, and depend on the lateral 
size of the membrane L: in the uj — > limit we find that 



T, 



off .direct — Tb OC g m f n wnile i e ff,curv — J-b oc q mm , 

5min = 2tt/L (Fig. [3] inset). By comparing the calculated 
and observed [l5| frequency dependence of T e g and the 
power spectral density [11], we can estimate the proper- 
ties of the active "motors" in this system: At ~ 100msec, 
Ponfo — (R/r) 2 , where r ~ lOOnm, R ~ 90ksT, and the 
motor density n = 1/r 2 . These parameters agree with 
the physical interpretation of the active force as arising 
from "pinching" of the membrane by a cytoskeleton net- 
work of spectrin filaments [3] • 

We demonstrate in Fig. [3] that a future measurement 
of T er T can be used to distinguish between different mod- 
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els of active forces in the membrane, and can therefore 
add important information about the physical nature of 
the motor producing these forces. In particular, a non- 
monotonic behavior will favor the curvature mechanism, 
while a simpler step-like behavior will support the direct 
force. For the curvature force, the effective temperature 
decreases with decreasing oj, and even approaches the 
equilibrium value at w — > for a large membrane domain 
(L — > oo). Another indication for an increase of effective 
temperature with frequency was found for a driven gran- 
ular system 17J. For the RBC membrane this behavior 



is driven by the fact that the curvature force couples 
to the fluctuation modes of the membrane through a q 2 
term 111 ] . which represents in g-space the force due to 
a localized induced curvature. This force therefore di- 
minishes in its relative amplitude as the wavelength in- 
creases, leading to the vanishing of the active component 
in the uj — > limit. We expect this feature to appear in 
many spatially extended systems where the driving force 
decreases with increasing wavelength. 

Conclusion. We presented a simple model for active 
systems, for which we can derive two measures to char- 
acterize its non-equilibrium nature. These two measures 
do not always agree. In biological systems the activity 
is often driven by multiple molecular motors that couple 
to internal degrees of freedom. In such systems we ex- 
plored the characteristics of the non-equilibrium fluctu- 
ations in the presence of multiple motors. We find that 
in the limit of many motors the kurtosis can return to 
the value of a Gaussian distribution, while the effective 
temperature may still exhibit strong frequency depen- 
dence. We showed that Gaussian distributions may arise 
for active systems even for a small number of motors, 
while large deviations from Gaussian distributions can 
be maintained even for large number of motors. Note 
that a non-Gaussian distribution is not a proof of non- 
equilibrium, as it could also arise due to nonlinearities in 
a mechanical system, such as position dependent damp- 
ing. The effective temperature and the kurtosis that we 
calculated are explicitly dependent on the number of mo- 
tors (N m ) and their intrinsic properties (/o,t, At). Our 
present analysis gives a detailed and general treatment, 
for any type of pulse-length distribution. 

Finally, we compared the results of our model with 
recent observations of ATP-drivcn activity in RBC, and 
demonstrated how they can give insight to the underlying 
active mechanism. In particular, we showed how funda- 
mental physical properties of the elusive molecular motor 
of the RBC membrane may be unraveled by comparing 
these observables with our theoretical model. Future ex- 
periments could use the calculated properties to better 
characterize the nature of the active forces in various cel- 
lular membranes. We expect our results to be useful for 
the analysis of other active systems, both biological @ 
and non-biological fl8L [l9j . From our model we reach 
the following more general conclusion: when a spatially 



extended system is driven by external forces, the effec- 
tive temperature defined through the FD relation can be 
a non-monotonic function of the frequency. If the cou- 
pling of the external active forces is stronger for smaller 
wavelengths, then the effective temperature may develop 
a non-monotonous dependence on frequency. 
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